HUPD-9907 

RESCEU-16/99 

UTAP-328/99 



The cosmological light-cone effect on the power spectrum 
of galaxies and quasars in wide-field redshift surveys 



Kazuhiro Yamamoto, Hiroaki Nishioka 
Department of Physics, Hiroshima University, Higashi- Hiroshima 739-8526, Japan. 

and 
Yasushi Suto 

Department of Physics and Research Center for the Early Universe (RESCEU) 
School of Science, University of Tokyo, Tokyo 113-0033, Japan. 



ABSTRACT 

We examine observational consequences of the cosmological light-cone effect on the 
power spectrum of the distribution of galaxies and quasars from upcoming redshift 
surveys. First we derive an expression for the power spectrum of cosmological objects 
in real space on a light cone, P^^Jyk), which is exact in linear theory of density 
perturbations. Next we incorporate corrections for the nonlinear density evolution 
and redshift-space distortion in the formula in a phenomenological manner which is 
consistent with recent numerical simulations. On the basis of this formula, we predict 
the power spectrum of galaxies and quasars on the light cone for future redshift surveys 
taking account of the selection function properly. We demonstrate that this formula 
provides a reliable and useful method to compute the power spectrum on the light 
cone given an evolution model of bias. 



Subject headings: cosmology: theory - dark matter - large-scale structure of the 
universe - galaxies: distances and redshifts - quasars: general 
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1. 



Introduction 



The importance of the cosmological hght-cone effect in the future redshift surveys, including 
the Two-degree Field (2dF) and the Sloan Digital Sky Survey (SDSS), has been well recognized 
recently (e.g., Matarrese et al. 1997; Matsubara, Suto, & Szapudi 1997; Nakamura, Matsubara, 
& Suto 1998; Moscardini et al. 1998; de Laix & Starkman 1998). In a series of previous papers 
(Yamamoto & Suto 1999, Paper I; Nishioka & Yamamoto 1999a, Paper II; Suto et al. 1999, Paper 
III), we have extensively explored and discussed various aspects of the cosmological light-cone 
effect on the two-point correlation functions, ■^^^(r), of objects at high redshifts. Extending the 
procedure described in these papers, we develop here a formula to predict the power spectrum of 
cosmological objects located on the light cone, P^^(fc). 

Inhomogeneities of the matter distribution might affect the light propagation, i.e., the 
gravitational lensing effect. In fact, the angular correlation functions are significantly affected by 
the gravitational lensing through the magnification bias. The spatial correlations which we study 
in the present paper, however, are fairly insensitive to the effect (e.g., Moessner, Jain, & Villumsen 
1998). Therefore we ignore lensing effect below. 

The plan of the paper is as follows. In §2, we define an estimator for the power spectrum 
constructed from a sample of cosmological sources on a light cone, which generalizes the 
conventional power spectrum defined on the constant-time hypersurface. By computing the 
ensemble average of the estimator, we derive the power spectrum on the light cone. A mathematical 
derivation is outlined in Appendix A. While this rigorous procedure assumes linear theory of 
density perturbations in real space, the effects of nonlinear gravitational growth and redshift-space 
distortion are important in practice (Papers II and III). Thus we incorporate those effects as 
described in §3. Then we apply our theoretical results to galaxy and quasar samples with the 
SDSS spectroscopic samples specifically in mind (§4). Finally §5 is devoted to discussion and 
conclusions. Throughout this paper we use the units in which the light velocity c is unity. 



In this section we describe our theoretical formulation for the power spectrum on a light-cone 
hypersurface which proceeds fairly parallel to Paper I. For simplicity, we focus on the spatially-flat 
Friedmann-Lemaitre universe, whose line element is expressed in terms of the conformal time r] as 



2. Power spectrum on the light cone: linear theory in real space 



2.1. Basic procedure 




We normalize the scale factor to be unity at present, 0(770) = 1- All cosmological objects observed 
from redshift surveys are located on the light-cone hypersurface of the space-time (|^), defined by 
an observer. We locate the fiducial observer at the origin of the coordinates {r] = r]Q, x = 0)) &iid 
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therefore an object at x ^'^^ V on the hght-cone hypersurface of the observer satisfies a simple 
relation r] = ijq — x- Then the (real-space) position of the source on the light-cone hypersurface 
is specified by (X)7)) where 7 is a unit vector along the line-of-sight. In order to avoid confusion, 
we introduce a radial coordinate r instead of X; and write the metric of the three-dimensional real 
space as 

dslc = dr^ + r^dnf2y (2) 

Let us denote the number density field of the sources on the light cone by n^*^ (r, 7) , which is 
simply related to the comoving number density of objects at a conformal time rj and at a position 

(X,7) as 

n^'^ir, 7) = n{r], X, 7) \rj^rjo-r, x-r ■ (3) 

Introducing the mean observed (comoving) number density noif]) at time r] and the density 
fluctuation of luminous objects A(?7,x,7), we write 

n(,V,X,7) = no{r])[l + A{ri,x,7)]- (4) 
Then equation (|3|) is rewritten as 

nLC(r, 7) = no(?7) [1 + A{7], X, 7)] \r,^r,o-r, x^r- (5) 

In what follows, we define 

n]^^{r) = no(r?) |^_^^,_, (6) 

for convenience. Note that nQ{r]) is different from the mean number density of the objects n{ri) at 
rj by the selection function which depends on the luminosity function of the objects and thus on 
the magnitude-limit of a survey (see Paper I and §4 below). 

Following the conventional treatment of the power spectrum in cosmology (e.g., Feldman, 
Kaiser, & Peacock 1994), we introduce the number density field of cosmological objects on the 
light cone: 

n^'^{r, 7) — n^*-'(r, 7) 



F(r,7) 



d-'rnl^ir) 



(7) 



where n^^(r, 7) denotes the number density of the synthetic catalog without structure, which has 
the mean number density nQ*-'(r). To be specific, it satisfies 

(nLC(ri,7i)nLC(r2,72)) = (nLC(ri, 7i)n^C(r2, 72)) = 4^(ri)4C(r2), (8) 

where we neglect the Poisson noise term. 

A power spectrum on the light cone, which can be estimated from a given survey, is written as 

/ d'k\T{k)\' 

P(/t)°b« = ^ (9) 

d^k 

Vs 
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where Vg is the volume of a thin shell with the radius k in the Fourier space, and ^(k) is the 
Fourier transform of equation (|^) : 



^(k) 



LC/ 



ikr 



/ 



^ 1 

2 



(10) 



Finally the power spectrum on the light cone, -PrhhC^), can be computed by taking the 
ensemble average of the above estimator (P): 



Pk%n{k) = = (|^(k)|2) 

d'n [ rfV2([nLC(ri,7i) -nLC(ri,7i)][nLC(r2,72) -nLC(r2,72)])e^'^-('-^-'-^) 



(11) 



For definiteness, we use a superscript LC explicitly to indicate the power spectrum of objects on 
the light cone throughout the paper. The power spectrum without the superscript denotes that of 
mass defined on the constant-time hypersurface. Also subscripts R and S indicate those in real 
and redshift spaces, and subscripts lin and nl refer to those in linear and nonlinear models. 

Assuming linear bias between density fluctuations of mass and luminous objects and also 



linear theory for growth of density fluctuations, we find that equation (11) reduces to the following 
expression: 



-PR,Un(^) 



— I 

2TTk Jo 



dnri J c?r2r2nQ^(ri)nQ^(r2)Di(r/o - ri)Di{r]o - r2) 
dkikiPuiiaiki,7]o)b{ki;rio - ri)b{ki;rjo - r2)IL{ri,r2,ki,k) , (12) 



where D\{rj) is the linear growth rate normalized unity at present, Di{r]Q) = 1, b{k;r]) is the linear 
bias factor, and 



U{ri,r2,ki,k) = 
with the cosine integral function Ci{x) defined as 



Ciiln - r2\\ki - k\) + Ciiln - r2\{ki + k)) 

+Ci{{n + r2)\ki -k\)- Q((ri + r2){ki + k)) (13) 



aix) 



°° , cost 

dt . 

t 



(14) 



The derivation of equation (12) from equation ( |l 1|) is quite tedious but straightforward which we 
outline in Appendix A. 
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2.2. Approximation to the exact formula 



While equation (12) for power spectrum on the light cone is exact as long as linear theory of 
density fluctuations and a linear bias model are adopted, it looks fairly complicated. Therefore 
it is instructive to consider a practical approximation to equation One can show that the 

function n(ri, r2, ki, k) (eq.||l^) is peaked around ri ~ r2 and /c ~ fci. In fact, if one replaces the 
function as 



U{ri,r2,ki,k) ~ 2iT5{ki - k)5{ri - r2) 



equation (|12|) reduces to 



^R,un(fc) - aR,iin(A:)PR,iin(A:,z = 0), 



with 



aR,lm(^) 



-1 



drr'^nY^{rfDi{r]Q - rYb{k;r]Q - rf 



(15) 



(16) 



(17) 



In an extreme case of no evolution without biasing, i.e., b{k;r]) = Di{r]) = 1, aYi,hnik) becomes 
unity, and fRnn(^) equivalent to i'(/c)R,iin as expected. Thus if the approximation (p^ is 
justified, ajifi^ik) quantifies the degree of the light-cone effect on a specified survey, which 
sensitively depends on the observed selection function and the assumed time-evolution of linear 
bias. 



Figure || compares the power spectra on a light cone, equations (12) and (16). It should 
be noted that P^^iy^ik) in the exact formula becomes significantly larger than PRjm{k) for 
k < l/r( 

^max) with r(zinax) being the depth of the survey volume. This artifact is not particular 
to the light-cone effect, but simply originates from the fact that the power on scales comparable 
to, or larger than, the finite size of the survey volume cannot be properly evaluated. 

Thus except for the range k <^ l/?'(-Zmax)) where the estimate of the power is fairly reliable 
for a given observational sample. Figure || shows that the approximation ([l6| ) reproduces the exact 
formula very accurately. Therefore in what follows, we will incorporate several important effects 
(redshift-space distortion, nonlinear density and velocity evolution, and the observational selection 
function) on the basis of the approximation. 

While the results in Figure || do not include the selection function, i.e., no(z)=constant. is 
assumed, we also examined the case with the proper selection functions for galaxies and quasars 
(§4), and made sure that our approximation agrees with the exact formula within several percents 
at worst. This level of accuracy is sufficiently good in the light of the other approximations 
employed in describing the nonlinearity of density and velocity fields, most notably, the uncertainty 
of the model of bias. 
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3. Redshift-space distortion and nonlinear evolution of density and velocity fields 

The results presented in the previous section are unrealistically simplified since we have 
neglected several important effects including the redshift-space distortion due to the peculiar 
velocity effect (Kaiser 1987), nonlinear evolution of mass density fluctuations and of peculiar 
velocity dispersions (finger-of-God) , evolution of the bias parameter (e.g., Fry 1996), and the 
observational selection function which depends on a specific survey strategy. In this section 
we describe our modeling to the first three effects which are the necessary ingredients in 
completing the theoretical predictions. The selection function is discussed in the next section. 



Strictly speaking, inclusion of these effects invalidates the exact derivation of the formula (12) 
unfortunately. Therefore we implement a phenomenological correction to these effects as described 
below, but we believe that this procedure is important and useful in practice, especially given the 
uncertainty of the bias evolution itself. 



3.1. Linear redshift-space distortion 



The cosmological observation is possible only in redshift space. The redshift-space distortion 
in linear theory, first discussed by Kaiser (1987), is easily included in our formula for power 
spectrum on the light cone in real space. In linear theory, the direction-averaged power spectrum 
in redshift space, Ps,hn{k,z) is related to that in real space, PR,hn{k, z), as 



Ps,im{k,z) 



l + ^P{z) + -^P'{z) 



R,lin 



{k,z). 



(18) 



The above relation (Kaiser 1987) assumes a distant-observer approximation, and scale-independent 
linear bias, b{z), as well. Then P{z) is defined by 



P{z) 



1 dlnDi 
b{z) din a 



(19) 



See Matsubara & Suto (1996), Hamilton (1998) and Paper III for extensive discussion on the 
redshift-space distortion. 

Substituting equation ([l8|) in equation ( [T^ ) yields the expression for P^](^{k) which is 
approximated as 



^s:nn(^) - as,iin(A:)-PR,iin(A:,z = 0), 
as,iin(fe) 



J drr^nQ^{rfDi{rjo - r)%{k; r]o - r)^ 
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J 


f drr^nY^{rf 



(20) 



(211 



A similar expression for the linear redshift-space distortion on two-point correlation functions is 
derived in Paper II. 
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3.2. Nonlinear evolution of density 

Nonlinear evolution of mass density field introduces the additional /c-dependence on the 
light-cone power spectrum. While the full description of the nonlinear gravitational evolution 
of density fields is almost impossible, fairly accurate fitting formulae for PR ni(^)-2) have been 
obtained (e.g., Peacock & Dodds 1994,1996). We attempt to correct for the nonlinear density 
evolution simply by adopting the fitting formulae, and then obtain an approximate expression for 
^R,nl(^) as follows: 

^R,nl(^) - aR,nl(fc)PR,lin(fc,^ = 0), (22) 



aR,ni{k) = 7 • (23) 



3.3. Nonlinear redshift-space distortion; finger-of-God effect 

Finally an effect of the nonlinear velocity field, finger-of-God, should be taken into account in 
order to make a testable theoretical prediction. Again this has been extensively discussed for the 
power spectrum on the constant-time hypersurface (Cole, Fisher & Weinberg 1994; Mo, Jing, &: 
Borner 1997; Paper III; Magira, Jing & Suto 1999). In particular. Cole et al. (1994) proposed the 
following correction for the finger-of-God effect: 



-Ps,nl(^,^) 
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i^R,ni(A;,^). (24) 



If the velocity distribution function is approximated by an exponential model with scale- 
independent one-dimensional dispersion of the peculiar velocity (7v{z), which is indicated 
observationally (Davis & Peebles 1983) and also from numerical simulations (Efstathiou et al. 
1988; Ueda, Itoh & Suto 1993; Magira et al. 1999), the above functions are explicitly given by 

_ aretm(KA/2) !_ 
'''''' " V2.-. +2 + ^2- 

C(„,) = ^(b(«)-^), (27) 

with k{z) = kav{z) / Hq. On large scales, (Ty{z) can be well approximated by a fitting formula 
proposed by Mo, Jing & Borner (1997). The validity and limitation of the approximation are 
discussed in detail by Magira et al. (1999). Combining equation ( p^ with those fitting formulae, 
one can compute the nonlinear power spectrum on the light cone in redshift space as 

P^%{k) ~ as,ni(fe)i^R,iin(A:,^ = 0), (28) 
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as,ni(^) 



A{k) + |/3B(k) + llf-C(K) 

6 5 



PRfin{k,z = 0) I drr^n^^ir) 



(29) 



4. Predictions of power spectra on the light cone for the future redshift surveys 

4.1. Selection function for galaxies and quasars 

In properly predicting the power spectra on the Ught cone, the selection function should be 
specified. In this subsection, we describe the selection functions of galaxies and quasars with the 
upcoming SDSS spectroscopic samples in mind. 

For galaxies, we adopt a B-band luminosity function of the APM galaxies (Loveday et al. 
1992) fitted to the Schechter function: 

= ^* A) " exp (- ( A ) , (30) 

with (j)* = 1.40 X IQ-^/i^Mpc"^, a = -0.97, and = -19.50 + Slogig/i- Then the comoving 
number density of galaxies at z which are brighter than the limiting magnitude -Bum is given by 

POD 

no(z,<Aim)=/ 0(L)dL = 0*r[(a + l,x(Aim,^)], (31) 

where 

2 



x{Bii^,z) = 



L* 

and r[z^, x] is the incomplete Gamma function. 



dL{z) 



1/i-iMpc 



102.2-0.4B,,.^ (32) 



For quasars, we compute a selection function for B-band magnitude limited samples following 
Paper I on the basis of the luminosity function by Wallington & Narayan (1993; see also Nakamura 
& Suto 1997). 



4.2. Models and predictions 

For definiteness, we consider SCDM (standard cold dark matter) and LCDM (Lambda 
cold dark matter) models, which have {Qo,Qa, h, as) = (1.0,0.0,0.5,0.6) and (0.3,0.7,0.7,1.0), 
respectively. These sets of cosmological parameters are chosen so as to reproduce the observed 
cluster abundance (Kitayama & Suto 1997). We use the fitting formulae of Peacock &: Dodds 
(1996) and Mo, Jing, & Borner (1997) for the nonlinear power spectrum -Pa,ni(fc) and the peculiar 
velocity dispersions cr„, respectively. 

Figure |2| plots the theoretical predictions for power spectra on the light cone. We adopt the 
B-band limiting magnitude 19 and 20 for galaxies and quasars, respectively, so as to roughly match 
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the SDSS and 2dF spectroscopic samples. For illustrative purposes, the results are divided by the 
corresponding linear power spectrum on the constant-time hypersurface in real space at z = 0, 
PR,im{k,z = 0). For simplicity we adopt a scale-independent linear bias model of Fry (1996); 
b{k,z = 0) = 1 and 1.5 for galaxies and quasars (p = 1 in eq.||3^ below). Purely linear theory 
predictions, ^'s'un(^)' plotted in dotted lines, which use the CDM transfer function by Bardeen 
et al.(1986; BBKS) and the linear redshift-space distortion ( |20| ) by Kaiser (1987). We show two 
different versions of P^'^iik) ; both adopt the Peacock & Dodds (1996; PD) nonlinear spectrum, 
one, plotted in dashed lines, uses linear redshift-space distortion, and the other (solid lines) uses 
the nonlinear model ([2^ ) with the velocity dispersion formula by Mo, Jing & Borner (1997). 

In linear regime, the light-cone effect suppresses the amplitude of spectrum in real space 
due to the smaller fluctuation at larger z unless the possible evolution of bias exceeds the linear 
growth rate of the mass fluctuations. The redshift-space distortion, on the other hand, tends to 
increase the observable power spectrum in redshift space. Thus the overall effect of the linear 
redshift-space distortion on the light cone is to increase (decrease) the amplitude of P^'i^ni^) 
galaxy (quasar) samples in a /c-independent manner. 

The situation becomes complicated in nonlinear regimes, which produces the additional 
fe-dependent behavior. While the nonlinear evolution of density field enhances the amplitude of 
both -PRni(^) -^s'ni(^) relative to their counterparts in linear theory , the suppression due 
to the finger-of-God is much stronger. Consequently P^nii^) smaller than -Ps'iin(^) ^' 
and in fact the suppression is stronger for larger k. Since recent numerical simulations (Mo, Jing, 
& Borner 1997; Jing 1998; Suto et al. 1999; Magira et al. 1999) confirmed that the nonlinear 
fitting formulae of density and velocity fields are accurate and reliable within several percents, 
these conclusions are generic as long as the bias is time-independent. The bias, however, is 
predicted to show strong time evolution. Furthermore the realistic bias is expected to be neither 
time-independent nor linear (Fry 1996; Mo & White 1996); it may even not be deterministic 
(Tegmark & Peebles 1998; Dekel & Lahav 1999; Taruya, Koyama &: Soda 1999). We will discuss 
this problem in detail elsewhere, and focus on the effect of time-evolution of bias in the next 
subsection. 



4.3. Dependence on time-evolution of bias 

As discussed in previous subsection, the most uncertain aspect of the theoretical predictions 
of P^'^{k) is the model of bias. Since it is unlikely that any reliable theoretical model for bias is 
established in near future, all that one can try is to explore the possible effects by adopting a 
simple parametric model. While several phenomenological models of evolution of bias have been 
proposed in the literature (e.g., Moscardini, et al. 1998; Matarrese, et al. 1997), we adopt the 
following simple model: 

biv) = l + j^^[bo-l], (33) 
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with 60 = b(rjQ) is the bias parameter at present. When the constant index p is unity, the model 
reproduces the time-evolution of the lowest-order bias coefficient which is discussed by Fry (1996) 
on the basis of the continuity equation of mass and galaxy density fields. 

Figure |^ plots aR^un/^g and as,iin/^o for p = 1 (Fry's model) and 2 (regarded as an extreme 
example for comparison) as a function of 60 • The selection functions for the galaxy and quasar 
samples are identical to those adopted in Figure |2|. For shallow samples like the SDSS galaxies, 
theoretical predictions are fairly insensitive to the applied model of evolution of bias; the 
cosmological light-cone effect on these samples is dominated by the redshift-space distortion. For 
deeper samples, however, the amplitude is very sensitive to the evolution of bias; especially when 
the bias evolves faster than the linear growth rate {p > 1), the redshift-space distortion becomes 
relatively negligible on linear scales. 

5. Conclusions 

In this paper we have developed a theoretical formulation to compute the power spectrum 
of cosmological objects on the light-cone hypersurface. On the basis of the exact formula which 
works in linear theory and in real space, we have obtained a useful approximate expression 
valid on scales less than the survey volume size (Further discussions for the validity of the 
approximation will be given elsewhere; Nishioka & Yamamoto 1999b). In linear theory, the 
light-cone effect simply suppresses (in general) the amplitude of the observable power spectrum 
(as in the case of the two-point correlation functions; see Papers I, II and III). Then we improved 
the approximate formula by including the nonlinear evolution of density fields, and linear and 
nonlinear redshift-space distortion phenomenologically but in a manner fully consistent with the 
recent numerical simulations. These nonlinearities produce the additional scale-dependence in the 
power spectrum on the light cone compared with those defined on the constant-time hypersurface. 

Applying this expression to the galaxy and quasar samples selected so as to match the 
upcoming SDSS spectroscopic samples, for instance, we have quantitatively evaluated the degree of 
the light-cone effect. With these example, we have demonstrated that we have developed a fairly 
reliable and useful method to compute the power spectrum on the light cone given an evolution 
model of bias. In fact, it is clear that the evolution of bias sensitively changes the behavior of 
the power spectrum on the light cone, especially for quasar samples. While our present analysis 
assumes the simplest bias model, linear bias, we plan to extend our formulation for the general 
bias model including the stochastic bias models (Dekel & Lahav 1999; Taruya, Koyama, & Soda 
1999) 

We thank Y.Kojima, T.Matsubara and M.Sasaki for useful discussions and comments, and 
Y.P.Jing and H.Magira for providing numerical routines to compute fitting formulae for nonlinear 
density and velocity fields. This research was supported in part by the Grants-in-Aid by the 
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APPENDIX 



A. Linear theory derivation of the exact formula of the Ught-cone power spectrum 

-^RUn(^) real space 

In this appendix, we outline the derivation of equation (^) from equation (|ll|): 

/ d'ri f d3r2([nLC(ri,7i)-nLC(ri,7i)][nLC(r-2,72)-nLC(r-2,72)])e*k.(r,-r.) 



With equations (|5[) and (^),the numerator of the integrand in equation (Kl) becomes 

([nLC(ri,7i) -nLC(ri,7i)][nLC(r2,72) -n^C(r2,72)]) 

= nJ^C (n ) nf^C (r2 ) ( A (j?o - r 1 , ri , 71 ) A (7?o - r-2 , r2 , 72 ) ) . ( A2) 

We expand the fluctuations A(?7,r, 7) in equation (^) as 

A(f?>X,7)=/ dky2^kimiv)ykim{x,l), (A3) 
Jo , — 

l,m 

where yklmix?!) is the normalized harmonics on the flat space: 

ykim{x,1) = \ -kji{kx)Yirain^), (A4) 
V vr 

with ji{x) being the spherical Bessel function and Yjml^-y) being the spherical harmonics on a unit 
two-sphere. To proceed further, we adopt two assumptions; a linear bias model in Fourier space: 

A,Uv)=bik;v)itiv), (A5) 
and linear growth of the density fluctuations: 

41(^) = 41(%)^i(^), (A6) 

with Di{r]) being the linear growth rate normalized to be Di{r]Q) = 1. In the above expressions, 
we denote the density fluctuations by with CDM density fluctuation specifically in mind, 

although not essential at all. 

Then the power spectrum of mass density fluctuations at present, -PR,iin(^) is defined 
through 

(A7) 

Substituting equations ( [A^ ) to ([A7|), one can write the third term in equation (|A^ ) as 
{A{r]o - ri,ri,7i)A(r/o - r2,r2,72)) 



2 

TT Jo 



dkik^PR,^iiri{h,rio)b{h;Vo - n)&(fc2;% - r2)Di{riQ - ri)Di(??o - ^'2) 



X Y.ji{k^n)3i{kiT2)YCr^{n^^)YU^^^) ■ (A8) 

Lm 



Using the expansion of the plane wave in terms of the spherical harmonics: 



„iU.r 



Im 



and equation (|A8D , equation (|Al|) is explicitly written as 



^Rain(^) 



n -1 



dkiklPRfi^{ki,r]o)b{ki;rjo - ri)b{ki;rio - r2)Di{r]o - ri)Di{rjo 



x^ji{hri)ji{kir2)Yi*^{n^,)Yim{^^2) 

l,m 

xA7tJ2 i'''jLAkn)YL,MA^^.)YL,M,{^i:) 

LiAh 

L2M2 



Integrating equation ( [AlOj ) over and il^j results in 



^R,lin(A:) 



-1 



drir'i I dr2ri4^(riX^(r2) 



vr Jo 



dkiklP^,^li^{ki,r]o)b{ki;7]o - ri)b{ki;r]o - r2)Di(r/o - ri)Di{r]o 



;(47r)2^jKfcin)jK^ir2)j/(A;ri)jKA:r2)l]m(^^k)^m(^^k)- 



Using the mathematical formula: 

I 



21 + 1 

m=—l 



Air 



one can simplify equation ( All| ) as follows: 



pLC 
-'R,lin 



drr'^n\'^{rf 



drir'l I dr2rlnY^{ri)n^^{r2)Di{r]o - ri)Di{'qQ 



LC/ 



vr JO 



dkiklPB^firiih,Vo)bih;r]o - ri)b{ki;r]o - r2) 



X ^(2/ + I)ji{kiri)ji{kri)ji{kir2)ji{kr2). 
I 



Applying the formula: 

1 1 1 

ji{aX)jiib\) = — / dxia"^ + 6^ - 2abx)-2 sin[A(a2 + b^ - 2abx)2]Pi{x) 
2\ J -I 



and then 



2n -I- 1 

5{X -X')=J2 -^Pn{x)Pn{x'), 
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equation ( [A13| ) is written as 



2 roo 



drirf J dr2r|no^(ri)no'^(r2)-Di(r?o - ri)Di{r]o - 

2 roo 

- / dkiklP^^xi^{ki,r]Q)b{ki]riQ-ri)b{ki]r]Q-r2)y^{2l + l) 
ttJo Y 

1 1 1 

X— — / dxi{rl + To — 2rir2Xi)^2 sm[ki{r'f + r2 — 2rir2Xi)^]Pl{xi) 
2ki J-i 

1 1 1 

^TTT / dx2irl + To - 2rir2X2)~2 sin[A;(r? + r| - 2rir2X2) 2]P;(x2) 
2fe 7-1 

driri y dr2rln^'^{ri)nY^{r2)Di{r]Q-ri)Di{rio-r2) 
dkikiP^,^\i^{ki,r]o)b{ki;r]o - ri)b{ki;r]o - r2) 



drr^n};'^{rf 



1 -1 



nk 

^ sin[/ci(rf + — 2rir2x)-^/^] sin[/c(rf + - 2rir22;)-'^/^] 



-1 



+ r2 — 2rir2X 



(A16) 



Finally changing the variable x to y = (rf + — 2rir2x)^''^, equation ( |A16 ) reduces to 



-fR,nn(^) 



drr'^nQ^{rf 



1 -1 



drm / dr2r2nJj'^(ri)no^(r2)i:>i(r/o - n)-Di(?7o - ^2) 



/"oo rv-i-\-T2 sill k'l^y sin 

X—/ dkikiP^ii,,{ki,7]o)b{ki;r]o - ri)b{ki;r]o - r2) dy . (A17) 

vrfc Jo "'Iri-ral 2/ 

Rewriting the last integral of the above expression in terms of the the cosine integral function: 

Ci{x) ^ - / dt^, (A18) 
equation ( [A17 ) yields the final expression ([T^), namely, 



-fR,nn(^) 



drr'^riQ^irf 

1 

2^ Jo 



1 -1 



drm J dr2r2n^'^{ri)n^'^{r2)Di{r]Q- ri)Di{riQ- r2) 
dkikiPK^liniki,T]o)b{ki;r]o - ri)b{ki; r]o - r2)Il{ri,r2,ki,k), (A19) 



where 



n(ri,r2, ki,k) = -Ciiln - r-2pi - k\) + a{\n - r2\iki + k)) 

+Ci{{n + r2)\ki -k\)- Ciiin + r2){ki + k)) . 



(A20) 
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Fig. 1. — Power spectra of cosmological objects on a light cone, Pj(^i^{k). Upper and lower pairs of 
curves correspond to the samples of survey depth extending up to z = 0.2 and 3.0, respectively (the 
exact formula eq.[12| in solid lines, and its approximation eq.||l^ in dashed lines). For comparison, 
the power spectrum of the mass fluctuations defined on the constant-time hypersurface at present 
is plotted in dotted line. For definiteness we adopt SCDM model in which Qq = 1, JIa = 0, h = 0.5 
and o"8 = 0.6 



-17- 



10' 



ISI 



CL- 




IO' 



N 



3 

Oh 



1 1 1 I I I I I I 1 1 1 I I I M 

Galaxy (B<19, z<0.2) \ 
SCDM " 



Quasar (B<20, z<5.0) 
SCDM 

b(k,z = 0)=1.5 



Fry's bias model 




1 1 1 I I I I I I 1 1 1 I I I I M 



Galaxy (B<19, z<0.2) 



7 



LCDM 

b(k,z)=1.0 

/ 



Pi':„,: PD 

P|c„i: PD+finger-of-God 

I I I 1 1 m l — I I I 1 1 MIL — I I I I M ill — I I I 1 1 11 1 



Quasar (B<20, z<5.0) 
LCDM 

b(k,z = 0) = 1.5 / 
Fry's bias model / 



J^Q-l I I I I I I I I I I 

10-2 10-1 . 10-' 



J I I 



k [h Mpc-'] 



k [h Mpc 



Fig. 2. — Power spectra on the light cone, P [k), for various models, divided by -PR,iin(^, z = 0), 
the corresponding linear power spectrum on the constant-time hypersurface in real space at z = 0. 
Upper panels assume a galaxy sample with B < 19 and z < 0.2, while lower panels are for 
a quasar sample with B < 20 and z < 5, both of which roughly correspond to the upcoming 
SDSS spectroscopic samples. Left and right panels plot the results in SCDM and LCDM models, 
respectively. Thick and thin lines represent P^'^{k) measured in redshift and real spaces. 
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Fig. 3. — aR,iin/^o CKS,iin/&o ^ function of the linear bias parameter at present bo. Solid and 
dashed lines represent the ratio measured in redshift and real spaces, respectively; p = 1 in thick 
and p = 2 in thin lines. The model with p = 1 corresponds to the bias evolution model of Fry 
(1996). 



